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ABSTRACT 

We have used the Atacama Pathfinder Experiment 12 m telescope to detect 
circumstellar CO emission from MP Mus (Kl IVe), a nearby (D ~ 100 pc), ac- 
tively accreting, ~ 7 Myr-old pre-main sequence (pre-MS) star. The CO emission 
line profile measured for MP Mus is indicative of an orbiting disk with radius 
~ 120 AU, assuming the central star mass is 1.2 M & and the disk inclination is 
i ~ 30°, and the inferred disk molecular gas mass is ~ 3M e . MP Mus thereby 
joins TW Hya and V4046 Sgr as the only late-type (low- mass), pre-MS star 
systems within ~100 pc of Earth that are known to retain orbiting, molecular 
disks. We also report the nondetection (with the Institut de Radio Astronomie 
Millimetrique 30 m telescope) of CO emission from another ten nearby (D ~ 100 
pc), dusty, young (age ~ 10 — 100 Myr) field stars of spectral type A-G. We 
discuss the implications of these results for the timescales for stellar and Jovian 
planet accretion from, and dissipation of, molecular disks around young stars. 

Subject headings: circumstellar matter — stars: emission-line — stars: individual 
(MP Mus) 



1. Introduction 

Over the past two decades, astronomers have identified a few hundred young (ages ~ 5 
to ~100 Myr) stars, including more than a half-dozen T Tauri and post-T Tauri associations, 
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within ~100 pc of Earth (IKastner et al.lll997l ; IZuckerman fc Song||2004j ; iTorres et al.l 12008 
and references therein). Most such stars have been identified on the basis of their char- 
acteristically large X-ray fluxes, thermal infrared emission from warm circumstellar dust, 
optical spectral features (p rincipally, strong Li abso rption lines), and/or distinctive space 
velocities (see discussion in IZuckerman fc Song 120041 ) . Based on the low frequency of mid- 
to far-infrared excesses among the nearby, young stellar groups, onl y a small percentag e of 
such stars evidently retain massive, dusty, circumstellar disks (e.g., iRebull et al.ll2008l ). in 
contrast to the large disk fracti ons of pre-ma in sequence (pre-MS) star populations associ- 



ated with molecular clouds (e.g.. IWyattll2008l ). This contrast reflects the relatively advanced 



ages of the known young st ars within ~ 100 p c and suggests that, while many may still be 
building terrestrial planets ( IMelis et al.ll2010l ). most have already moved beyond the epoch 
of active giant planet formation. 

Examples of actively accreting, lo w-mass pre-MS st ars (i.e., classical T Tauri stars; 
cTTS) within ~100 pc are rarer still ( ITorres et al.ll2008l ). The archetype of these nearby 



cTTS systems — which are characterized by their unusually large Ha emission equivalent 
widths — is the intensively studied TW Hya. At a distance of 56 pc, TW Hya is among 
the closest known cTT£0. With an estimated age of ~ 8 Myr ( Zuckerman fc Song 2004 ; 



Torres et al.ll2008l ). TW Hya also remains one of the oldest known examples of a cTTS. In 



V4046 Sgr (age ~ 12 Myr, distance 72 pc 
erences therein) and th e cTTS MP Mus 
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20071 ; ICortes et al.ll2009l ). Estimates of the age and dista nce of the latter star range, respec- 



tively , from 6 Myr to 17 Myr and from 86 pc to 103 pc ( jMamajek et al.ll2002l ; ITorres et al. 



20081 ). The younger age (hence larger distance) ap pears more accurate, given the Li absorp- 
tion li ne strength of MP Mus ( jWeise et al.l 120101 ) and its likely association with the e Cha 
group (ITorres et al.ll2008l ). High- resolution (gratings) X-ray spectra of TW Hya, V4046 Sgr, 
and MP Mus reveal evidence f or a significant emission contribution from accretion shocks 



as opposed to coronal act ivity (IKastner et al.l 120021 ; IStelzer fc Schmittl 120041 ; iGiinther et al. 



20061 ; lArgiroffi et al.l 120071 ) . This is consistent with observations demonstrating that the Ha 
and UV emission from each of th ese three stars is stronger than that typically associated 
with pure chromospheric activity ( IHuenemoerder et al.l 120071 ). 



Radio molecular line emission studies of circumstellar planet-forming disks provide 
unique insight into the Jovian planet, Kuiper Belt, and comet formation zones within 



2 The TW Hya Association (TWA) members Hen 3-600 (D ~ 4 5 pc; IHuenemoerder et al.l 120071 and 



references therein) and TWA 30AB (D — 42 pc; Looper et al. 20 1 Oa llh ) are the closest known examples of 
actively accreting, low-mass stars. 
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the outer regions (tens to hundreds of AU in radius) of p re-MS circumstellar disks (e.g 



Zuckerman et al.l Il995l ; iDutrey et al.l 119971 ; iThi et al.l 12004 ). Thus far, however, only four 



star-disk systems within ~ 100 pc of Earth have been detected in molecular emission lines 
with a radio telescope: the aforementioned ~ 10 Myr-old, K-type cTTS systems TW Hya 
(IZuckerman et al.lll995r iKastner et al.lll997l ) and V4046 Sgr (iKastner et al.ll2008bl). and the 
A-typ e stars 49 Cet and HD 141569 (ages ~ 20 and ~ 5 Myr, respectively; IZuckerman et al. 



19951 ). Detections of a suite of molecular transitions toward TW Hya established that i t 
possesses a rich molecular disk viewed nearly pole-on (IKastner et al.l Il997l ; IThi et al.l 12004 ) , 
while the double-p eaked molecular line profiles characteristic of Keplerian rotation observed 
toward V4046 Sgr (IKastner et al.ll2008bl ) similarly demonstrated that this close (period 2.4 d) 
binary system is orbited by a circumbinary molecular disk. These single-dish radio molec- 
ular line spectroscopy results have subsequent ly been confirmed v i a radio interferometri c 
(Submillimeter Array) imaging of CO emission (iQi et al.ll2004l . 120061 ; iRodriguez et al.ll2010l ). 



In an effort to similarly characterize the gas mass, dynamics, and chemistry within the 
disk of MP Mus (spectral type Kl IVe), we searched for, and detected, CO emission from 
MP Mus with the Atacama Pathfinder Experiment (APEX) 12 m In this paper, 

we report on these results and their interpretation. We also report the nondetection of CO 
emission from a sample of ten other nearby (D ~ 100 pc), dusty, young (ages ~ 10-100 Myr) 
stars of spectral type A-G with the Institut de Radio Astronomie Millimetrique (IRAM) 30 
m telescope. 



Observations and Results 



2.1. IRAM 30 m telescope CO observations of nearby, dusty young stars 

We observed the ten starfl in Table ODin the 230.538 GHz J = 2 -> 1 transition of 12 CO 
with the IRAM 30 m te lescope in August 2009. These stars (Table [T]) were chosen from the 
list of iRhee et al.l (120071 ) on the basis of large infrared excesses (indicative of disk dust masses 
M& ~ 0.05M e ) and estimated ages ~ 100 Myr. We observed all target stars using IRAM's 



3 APEX is a collaboration between the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 

4 During the 2009 August IRAM 30 m observing run we also obse rved the young late- A/early- F sta r 
HR 8799, which i s orbited by mu ltiple giant planets and debris disks ( Marois et al. 20081 : Su et al. 20091) . 
As mentioned by Su et al. ( 20091) . our 12 CO(2-l) spectrum established that bright interstellar emission is 
present at the radial velocity of the star, precluding any attempt to constrain its circumstellar gas mass via 
single-dish CO measurements. 
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EMIR receivers and VESPA autocorrelator, with velocity resolution of 0.8 km s~ x . Spectra 
were acquired in wobbler switching mode, resulting in flat baselines for individual spectra. 
Total on-source integration times were between 1 and 2 hours per star. The weather was 
generally good (r 2 25 ~ 0.25-0.4) throughout the observations, with time-averaged system 
temperatures in the range 220-510 K. Pointing and focus were checked (against standard 
pointing sources and planets) every ~2 hours. Typical pointing errors were ~ 3", i.e., ~l/4 
beamwidth (FWHP 12" at 230 GHz). 

We used the CLASS^ radio spectral line data reduction package to sum over individual 
spectral scans obtained for each star, and then to subtract a linear-fit baseline from each of 
the resulting integrated spectra, calculating channel-to-channel noise levels in the process. 
No CO emission lines were evident in any of the resulting spectra obtained for the ten sample 
stars. Upper limits on main-beam brightness temperature T m t> (assuming F e ff = 0.94 and 
Vmb = 0.52 for measurements at 230 GHz with the IRAM 30 m telescope) and velocity- 
integrated CO line intensity / are reported in Tabled] as the values 3<7t and 30"/, respectively, 
where ot is the efficiency-corrected rms channel-to-channel noise in the CO spectrum and 
a i = <Jtv AvSv, with velocity resolution 5v = 0.8 km s _1 and an assumed linewidth Av = 3 
km s _1 . 



2.2. APEX 12 m telescope CO observations of MP Mus 

Service-mode APEX 12 m telescope observations of MP Mus ( J2000 coordinates a =13:22:07.55, 
S =-69:38:12.2) in the 345.796 GHz J = 3 -»■ 2 transition of 12 CO were performed on 15 and 
16 April 2010 UT. The receiver was SHFI/APEX2, using the facility FFTS backend with 
velocity resolution 0.21 km s _1 . Spectra were acquired in position switching mode using the 
wobbler. Total on-source integration time was 138 min. The weather was generally good 
throughout the observations (T225 ~ 0.2; pwv 1.8 mm; time-averaged system temperature 
~ 500 K). The best focus was established at the beginning of each day's observations, and 
regular pointing checks (using nearby planets and quasars as references) indicated pointing 
errors were ~ 2", i.e., ~ § beamwidth (FWHP 17.5" at 345 GHz). APEX spectra of MP 
Mus were reduced using CLASS, as described in §2.1. 

The resulting, integrated, baseline- and beam-efficiency-corrected spectrum (assuming 
Vmb — 0.73 for measurements at 345 GHz with the APEX 12 m telescope) is presented in 
Fig. [TJ This spectrum clearly reveals the detection of 12 CO(3-2) emission from MP Mus. 
From a Gaussian fit to the line profile, we find the central velocity to be +3.9 ±0.3 km s _1 



5 See |http: / / iram.fr/IRAMFR/GILDAS71 
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with respect to the local standard of rest (LSR). This translates to a heliocentric velocity of 
11.3 ± 0.3 km s _1 , which is consistent with the systemic veloc ity of MP Mus as determined 



from optical spectroscopy (11.6±0.2 km s 1 ; lTorres et al.ll2006l ). The integrated line intensity 



obtained from the best-fit Gaussian is 0.21 ± 0.03 K km s 1 . 



3. Analysis 



3.1. MP Mus: constraints on molecular disk radius and mass 



Given the signal-to-noise ratio of the CO (3-2) data obtained for MP Mus (Fig. [TJ, we 
cannot unambiguously determine the underlying emission line profile. However, the steep 
sides and possible central valley (i.e., double-peaked appearance) of the MP Mus CO emission 



line a re features expected in the case of an orbiting molecular disk (e.g.. lBeckwith fc Sargent 



19931 ). Hence, we fit the profile with a parametric representa tion of the line pr ofile predicted 



by the Keplerian disk model of iBeckwith fc Sargentl (119931 ) as described in iKastner et al. 



(j2008al ). Provided the observed line profile is well fit, this method yields an estimate of outer 
disk rotation velocity — where Vd is equivalent to the half- value of the velocity separation 
of the red and blue peaks in the double-peaked line profile predicted for Keplerian rotation 
- as well as measures of the slopes of the inner and outer portions of the line profiles (pd 
and q, respectively). The value of q serves as an indication of the slope of the disk radial 
temperature profile (T oc r~ q ), while pd indicates the degree of central filling of the line 
profile by emission from gas at low radial velocities (e.g., for a nearly edge-on disk, pd = 1 
would correspond to a sharp outer edge and values pd < 1 would indicate lack of a sharp 
edge) . 

The observed line profile is overlaid with the best-fit model profile in Fig. [IJ The 
comparison confirms that the CO line emission profile of MP Mus is consistent with that 
of a molecular disk in Keplerian rotation. The best-fit Keplerian model parameter values 
(and formal uncertainties on these values) are Vd = 1.6 ± 0.16 km s _1 , q = 0.5 ± 0.2, and 
Pd = 0.25 ±0.1. The best-fit value of q is compatible with standard models of Keplerian 
disks 



(e.g. 



Beckwith k Sargentl 119931 ) . 



Given an estimate for the mass of MP Mus, M*, the best-fit value of Vd can be used 
to place constrain ts on the molecular dis k outer cutoff radius R ou t,co for an assumed disk 



inclination % (e.g.. IZuckerman et al.ll2008l . and references therein) . Adopting M ± = 1.2 M® 
(based on comparison with pre-MS evolutionary tracks; iMamajek et al.ll2002l ) and % = 32° 



(based on the morp hology of the disk as detected in coronagraphic imaging of starlight 



reflected off of dust; ICortes et al.ll2009l ). the result Vd = 1.6 km s 1 implies R, 



out,CO 



120 
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AU. This is somewhat sm aller than th e oute r disk radius inferred (also on the basis of 
coronagraphic imaging) by ICortes et al.l (120091 ) . i.e., R ou t,d ~ 170 AU, assuming a distance 
to MP Mus of 86 pc. However, given the uncertainties in our line profile analysis as well as in 
the mass, age, and distance of MP Mus (§4) — and the possibility that the (3-2) transition 
of CO is not well excited within the outermost regions of the disk — it appears our estimate 
of the radius of the molecular disk is in reasonable agreement with the dust disk radius as 
measured coronagraphically. 

The best-fit Keplerian model 12 CO (3-2) line profile area is slightly larger than that 
obtained from the Gaussian fit (§2.1), i.e., 0.26 ± 0.03 K km s _1 . This is a factor ~ 7 
weaker than the inte grated 12 CO (3-2) line intensity of TW Hya as measured with the JCMT 
( IKastner et al.lll997l ). suggesting that the intrinsic total 12 CO(3-2) line intensity of MP Mus 
is a factor ~ 3 smaller than that of TW Hya, after accounting for the difference in JCMT and 
APEX beam sizes and assuming MP Mus lies at ~ 100 pc (§4). The measured total APEX 
12 CO (3-2) line intensity of MP Mus is equivalent^ to an integrated line flux of 10.5 ± 1.5 
Jy km s~ l . This result can be u sed to estimate the molecular gas mass of the MP Mus 
disk (e.g., IZuckerman et al.ll2008l and references therein) — albeit with large uncertainties, 
especially given that thus far we have no measurements of other CO isotopes and (hence) 
no means to estimate the 12 CO(3-2) line optical depth (r 32 ) other than the foregoing rough 
comparison with TW Hya. Adopting the standard assumptions of a CO:H2 number ratio 
of 10~ 4 and molecular gas excitation temperature T ex = 40 K, and assumin g T32 = 3 (i.e., 
a factor ~ 3 smaller than that inferred for TW Hya; IKastner et al.l l2008bl and references 
therein) we obtain a molecular (H 2 ) gas mass of ~ 9 x 10~ 6 M Q (i.e., ~ 3M ffi ). This inferred 
disk gas mass, which scales approximately linearly with T32 if the 12 CO (3-2) is optically 
thi ck, is somewhat l ess th an the circumstellar dust mass of 17 M® determined for MP Mus 
by lCarpenter et al.l (j2005[ ) on the basis of the star's 1.2 mm flux. 



3.2. Other nearby, dusty, young stars: gas mass upper limits 

We used the same method just described for MP Mus to estimate molecular gas mass 
upper limits for the dusty disks of the ten stars observed with, but not detected by, the IRAM 
30 m. The results, obtained under the assumptions CO:H 2 = 10~ 4 , T ex = 40 K, and r 2 i = 1, 
are reported in Table [TJ These upper limits are of course subject to the same uncertainties 
just described in the case of MP Mus. Nevertheless, the results in Table [T] indicate that, for 
these ten stars, the disk gas masses are constrained to be of order of (or smaller than) a few 
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times the dust masses inferred from far-IR (IRAS) data, provided the CO emission is not 
optically thick or that gas-phase CO in these disks is not severely depleted relative to H 2 . 



Discussion 



With our CO detection of a disk orbiting MP Mus, this star joins ~8 Myr-old TW 
Hya and ~ 12 Myr-old V4046 Sgr as the only late-type (therefore low-mass) pre-MS stars 
within ~ 100 pc of Earth that are known to retain orbiting molecular disks. In all three 
cases, the circums tellar gas masses inferred from far-IR and radio emission lin es — a few to a 
few hundred M e (IKastner et al.lll997l : iRodriguez et al.ll2010i : iThi et al.ll2010l . and references 
therein) — are similar to the dust masses determined from far-IR and (sub)mm photometry. 
This implies dust-to- gas ratios near unity in the outer (~ 10 AU) regions of these disks 
and/or, in the case of gas mass estimates based on radio CO lines, that gas-phase CO is 
severely depleted relative to H 2 . 

The overall similarity of MP Mus to the T W Hya and V4046 Sgr star-disk systems favors 
the younger end of t he range of a ges es timated to date for MP Mus. This is consistent with 
the recent results of IWeise et al.l (120101 ). who estimated an age of 7 ± 3 Myr based on the Li 
absorption line streng th of MP Mus. An age of ~ 7 Myr would, in turn, imply the distance to 
MP Mus is ~ 100 pc ( Torres et al.l2008 ). For such a distance, the mass of MP Mus could be 
~ 20% larger than previously determine d (i.e., closer to 1.4 M ), depending on the adopted 
theoretical pre-MS evolutionary tracks (IMamajek et al.ll2002l . and references therein). If so, 
the CO and dust disk radii would be correspondingly larger than the estimates stated in 
53.1. 



Given the relatively advanced pre-MS evolutionary states of T W Hya, V4046 Sgr, and 
MP Mus, their circumstellar environments likely are in transition from predominantly pri- 
mordial material to "debris" disks that are populated by d ust grains generated via collisions 
between rapidly growing planetesimals (jCortes et al.ll2009l . and references therein). The stars 
for which we failed to detect molecular gas (Table []JJ likely fall into this latter category; all 
of the Table [1] stars have disk dust masses much smaller than those of TW Hya, V4046 Sgr, 
and MP Mus, and all have ages >10 Myr. The comparison of the residual gas masses of the 
three nearby cTTS with the gas mass upper limits in Table [1] thereby further reinforces the 
notion that the gaseous disks required to spawn Jovian planet s dissipate rapidly, i.e., within 
the first few Myr of pre-MS evolution (jZuckerman et al~ 1995 ). 



"Transition disks" usually display evidence for inner (radii ~ a few AU) disk clear ings 
that may be the results of recent or ongoing planet formation (e.g., ICalvet et al.l 120051 ). In 
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this respect, MP Mus, TW Hya, and V4046 Sgr are typical of transition disk systems; all 
three lack the significant near-IR excesses, indicative of hot dust in inner disks, tha t are 



characteristic of rapidly accreting cTTS in molecular clouds (e.g.. iMeyer et al.l 119971 ) . At 



the same time, like certain cTTS in Taurus that display transition disks (DM Aur and 



GM Aur; ICafvet et al.l 120051 ). all three nearby, late- type stars with detectable circumstellar 
CO emission also display spectral evidence for ongoing accretion — evidence that (to our 
knowledge) is lacking in the cases of the Table [1] stars. The fact that TW Hya, V4046 Sgr 
and MP Mus are evidently still accreting from their "transition d isks" may indicate that they 



are orbited by young planets (see discussion in lSalyk et al.ll2007l ). Regardless, the detections 



of CO in the dusty disks of MP Mus, TW Hya and V4046 Sgr suggest that the presence of 
residual disk molecular gas is closely linked to (likely enables) sustained accretion activity 
at advanced pre-MS ages. 

This research was supported via NASA Astrophysics Data Analysis grant NNX09AC96G 
to RIT and UCLA. The authors wish to thank the staff of the APEX telescope for their 
expertise in carrying out the observations of MP Mus. 
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MP Mus: C0(3-2), APEX 12m 



0.08 



0.06 



0.04 



0.02 



0.00 



-0.02 



J 



/r 
in 



Jl 



i i, 

i 



[fi 



ru 


-1- 




lr U 







-0.04 



■10 



10 
V sys [km/s] 



20 



Fig. 1. — 12 CO (3-2) line profile of MP Mus as measured by the APEX 12 m telescope and 
spectrometer (solid line) overlaid with best-fit Keplerian disk model profile (dashed line). 
Ordinate is velocity with respect to the local standard of rest and abscissa is main-beam 
brightness temperature. The spectrometer data have been Hanning smoothed to a velocity 
resolution of 0.4 km s -1 . The rms noise level of the spectrum is 13 mK per 0.4 km s -1 bin. 



Table 1: IRAM 30 m CO Observations of Nearby, Dusty Stars 
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-35:41:43.6 


A2 


111 


10 


0.2 


207 


1.1 


8.4 


< 4.8 


191089 


20:09:05.21 


-26:13:26.5 


F5 


54 


< 30 


0.03 


132 


0.69 


5.4 


< 0.7 


221853 


23:35:36.15 


08:22:57.4 


F0 


71 


< 100 


0.05 


57 


0.30 


2.4 


< 0.6 



a) Spectral types, distances, ages, and estimated dust masses as listed in lRhee et al.l ([20071 ) . 



